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Abstract 

The evolution of a quasi-isolated finite quantum system from a nonequilibrium 
I initial state is considered. The condition of quasi-isolation allows for the description 

of the system dynamics on the general basis, without specifying the system details 
and for arbitrary initial conditions. The influence of surrounding results in (at least 
partial) equilibration and decoherence. The resulting equilibrium state bears infor- 
mation on initial conditions and is characterized by a representative ensemble. It is 
shown that the system average information with time does not increase. The partial 
equilibration and non-increase of average information explain the irreversibility of 
time. 
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1 Introduction 



The feasibility of producing various kinds of finite quantum systems has sparked recently 
a strong interest in studying nonequilibrium properties of such objects. As examples of 
finite systems, it is possible to mention spin assemblies, quantum dots and wells, and 
trapped atoms. The current studies follow mainly two trends. The first direction, started 
yet by Schrodinger [1,2], considers thermalization of a non-isolated system coupled to a 
giant thermal reservoir. Related to this is the case of local equilibration, when a finite 
part of a large system is considered, with the rest of the system playing the role of a 
giant reservoir. The second track, pioneered by von Neumann [3,4], explores relaxation 
of isolated quantum systems to quasi-stationary states characterized by ergodic averages. 
The relaxation of composite finite systems also pertains to this track [5-7]. There have 
appeared a number of papers on these topics of equilibration. Extensive list of literature 
can be found in the recent review articles [8-13]. 

The aim of this paper is to study the third way, when a finite quantum system is 
quasi-isolated so that, though there exists some surrounding, but it is such that it does 
not disturb the system. In mathematical terms, this means that the matrix representation 
of all operators from the algebra of local observables is not infiuenced by the presence of 
the surrounding. It turns out that such quasi-isolated systems enjoy interesting specific 
properties, essentially differing them from both non-isolated as well as isolated systems. 



2 Quasi-isolated quantum systems 

Let us consider a system described by a Hamiltonian Ha acting on a Hilbert space T-La- 
And let there exist a surrounding characterized by a Hamiltonian Hb on a Hilbert space 
Hb- The total Hamiltonian is the sum 

Hab — Ha + Hb + Hint , (1) 

where H^nt is an interaction part. Hamiltonian (1) does not depend on time and acts on 
the space 

nAB = nA^HB. (2) 

The system statistical ensemble is the pair {T-Labi pAsii)} of space (2) and statistical 
operator 

pAB{t) = UAB{t)pABUXB{i) > (3) 

with the evolution operator 

Uab = exp{-iHABt) . (4) 

Here pab = Pab(O). The statistical operator is normalized, Ti ABPAB{t) = 1, where the 
trace is over space (2). 

The quantities of interest are the averages of the operators A from the algebra of 
local observables acting on the space Ha- Such averages correspond to the measurable 
quantities 

( A{t) ) = TvABPABit)! = TrAPAit)A , (5) 

where 

PA{t) = TYBPAB{t) . (6) 
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The system, by definition, is quasi- isolated from the environment when the system 
Hamiltonian is conserved, such that 

[Ha, Hab] = . (7) 

As a consequence, we have 

[Ha, Hb + Hin^ = [Ha, Hi^t] = . 

The quasi-isolation condition (7) guarantees that the matrix representation of all system 
operators from the algebra of local observables remains unchanged in the presence of 
the environment. In that sense, the latter does not disturb the system. Condition (7) 
slightly differs from the notion of a nondemolition measurement [14]. According to the 
latter, a nondemolition measurement of an observable, represented by the operator A, is 
such that this operator is an integral of motion, [A, Hab] = 0. Particular systems in an 
environment satisfying Eq. (7) have been considered, assuming that this environment is a 
bath composed of an infinite number of harmonic oscillators [15-19]. But here, we do not 
specify the environment. It will be shown below that the quasi-isolation condition (7), as 
such, is sufficient for deriving the general properties of equilibrating quantum systems. 
Let the eigenproblem for Ha be of the form 

HA\n)^E„\n) . (8) 

Then, because of commutator (7), the eigenproblem for Hamiltonian (1) reads as 

Hab \nk) ^{Er, + Snk) \nk) , (9) 
where \nk) = \n) l^\k) and Snk is the eigenvalue of the problem 

{HB + Hir,t)\nk) ^Enklnk) . (10) 

The quantities of interest are the averages of the operators A from the algebra of 
local observables acting on the space Ha- Such averages correspond to the measurable 
quantities (5) describing the system behavior. Using the basis composed of the vectors 
\nk) gives 

{A{t))^J2P^nmnm, (11) 

mn 

where A^^ = (m|74|n). The reduced density matrix 

Pinit) = { m \pA{t)\ n ) = Y,Pmnk{t) (12) 

k 

is the sum of the terms 

Pmnk{t) = { mk\ pAB{t) \ uk ) = pmnk{^) exp{-i(a;^n + £mnk)t} , (13) 

with the initial values 

Pmnifc(O) = {mk\ Pab(O) \nk) , (14) 
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and with the transition frequencies 
As follows from Eqs. (15), 

^nn — , ^nnk — . 

Hence the diagonal elements of matrix (12) do not depend on time: 

k 

Then, separating in average (11) the diagonal and nondiagonal terms yields 

( Mt) ) = ^Pnn(0)^„n + Yl pLit)Anm ■ 

The density matrix (12) can be written as 

Pmnit) = Rmn{t) exp{-iuJmnt) , 

where 

Rmn{t) = J^Pmnfc(0)exp(-i£„„fci) . 
k 

Introducing the density function 

gmni^) = y^Pmnfc(0)^(g " ^mnk) 
k 

makes it possible to rewrite Eq. (20) as the integral transformation 

/oo 
■oo 

The density function (21) is normalized so that 

/oo 
gmn{£) de = ^Pmnfe(O) = Pmn(O) ■ 
k 

This allows us to represent it as 

9mn{^) — Pmn(0)Pmn(^) ) 

with the distribution function normalized to one. 



/CO 
-oo 



Thence Eq. (20) can be written as 

Rmnit^ Pmn(0)-Omn(^) ) 



with the decoherence factor 



/oo 
Pmn{e)e-''' de . (27) 
-oo 

Consequently, matrix (19) becomes 

PLit) = Pmn{t)Dmn{t) , (28) 

where 

Pmn{t) = Pmn{0) exp{-iUmnt) ■ (29) 

As a result, for average (18), we have 

( A{t) ) = J2 Pnn{0)Ann + J] Pmn{t)AnmDmn{t) ■ (30) 

3 Possible types of equilibration 

If the surrounding is finite, so that the spectrum, defined by eigenproblem (10), is discrete, 
then the density function (21) is a sum of delta functions. Therefore factor (27) is a sum 
of exponentials and average (30) is a quasi-periodic function. In such a case, there is 
no equilibration in the strict sense, as far as any initial value of the observable will be 
reproduced after the Poincare recurrence time. However, if the number of degrees of 
freedom, characterizing the surrounding, is sufficiently large, the recurrence time can be 
very long [12]. Then the system tends to a quasi-equilibrium state, in which it lives during 
most of the time in the temporal interval before the recurrence time. 

But if the surrounding is so large that can be treated as infinite, then the multi-index 
k becomes continuous and the sums over k are to be replaced by integrals. In that case 
the density function (21) can be measurable, similarly to the properties of densities of 
states for large statistical systems. If so, the related state distribution, introduced in Eq. 
(24), is also measurable. Then the following general statement holds true. 

Theorem 1. // the state distribution Pmni^) is measurable, then the quasi-isolated 
quantum system equilibrates in the strict sense, so that 

lim(i(t)) = Vp„„(0)A„, . (31) 

t— >-oo ' 
n 

Proof. The state distribution, by assumption, is measurable and, by definition (25), 
it is integrable. Then, according to the Riemann-Lebesgue lemma [20], for integral 
(27), one has 

lim Dmnit) = . 

>oo 

The matrix Pmn{t) depends on time through the bounded function exp(— icu^^t). Since, 
by the Riemann-Lebesgue lemma, the decoherence factor tends to zero, we come to the 
limiting equality (31). 
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As an illustration, showing how the density function (21) can become measurable, let 
us treat the index k as continuous, so that function (21) be represented by the integral 



gmn{£) = J PmnkiO)S{e - Emnk) dll{k) , 



with a d-dimensional differential measure 

27r'^/2 

= wm * * • 

corresponding to a spherically symmetric surrounding. 

The delta-function of a function f{x) is defined by the expression 



5{x- Xj) 



in which Xi are simple zeros given by the equation f{xi) — 0, with the derivatives f'{xi) 
being nonzero. Then, for our case, we get 

r/ , N ^{k - ki) 

0{e — Emnk) — / ^ 



mnki 



where k^ — ki^ni^) is defined by the equations 

^ _ ^ — n ir' = ^^rnnkj I ^ 

'^mnki — ^ ) '^mnki — 7" ^ ■ 

This reduces the density function to the measurable form 

(0) ^ 



mnki I 

The second term in Eq. (30) is caused by quantum coherence. Therefore, its disap- 
pearance because of the action of environment, as in limit (31), is named the environment 
induced dccoherencc [21-23]. 

The concrete expression for the decoherence factor (27) is defined by the state dis- 
tribution Pmn{^) induced by surrounding. Keeping in mind a large number of random 
perturbations, produced by the surrounding, the induced distribution, in view of the cen- 
tral limit theorem, can be modeled by the Gaussian form, with a standard deviation 

^mm ^ 

Pmn(£) = CXp 



As a result, the decoherence factor is also Gaussian 



D^nit) = exp ^ - i {a^^tf \ . (32) 



In the case of the Lorentz distribution 



Prani^) 



+ lln) ' 



we have the exponential form 



D^Jt) = eM-lmnt) . (33) 
If the induced density is of the Poisson type 



Pn.nl^J = ^ expl -- 



we get the decoherence factor 



When the surrounding is modeled by a uniform distribution 

= 7^ e>(Amn - £)e(A^„ + e) , 

the decoherence factor is of power law with oscillations: 

( ^mn^) 



^L(^) = : 7 ^ ■ (35) 



Under different circumstances, one can meet different behavior of the decoherence 
factor. For instance, the relaxation to an equilibrium state can be Gaussian [24], or ex- 
ponential [15,19,25], or of power law [16]. Generally, the influence of surrounding can 
be represented by a linear combination of the above types of induced densities. Hence, 
the decoherence factor can also consist of several parts. For example, it can be a com- 
bination of the Gaussian and exponential terms [11,17,18]. In any case, if the density 
Pmn{^) is measurable, then a quasi-isolated quantum system equilibrates to a stationary 
state characterized by the diagonal elements p„ = p„„(0), whose deflnition involves in- 
formation on the initial state of the system, which needs to be given. The most general 
method of characterizing a stationary statistical system is by constructing the corre- 
sponding representative ensemble by maximizing the system entropy, under prescribed 
additional constraints [26-28], which is equivalent to the minimization of the information 
functional [29,30]. One obvious constraint is the normalization of p„. Other constraints 
are given by the average values of constraint operators Ci, usually corresponding to some 
observable quantities Cj and initial conditions [31]. 

A quasi-isolated quantum system can also thermalize, when on the manifold of initial 
conditions there exists a basin of attraction, such that limit (31) docs not depend on the 
particular choice of initial conditions from this basin of attraction [12]. A simple case is 
when the system initial state is an eigenstate, say so that Pnn(O) = Snj, then limit (31) 
is exactly Ajj. Suppose now that the initial state is prepared such that Pnn(O) is nonzero 
only for a window of the indices n e Nj around j, where the eigenvalue of Ha in Eq. 
(8) is Ej. And let this window be so narrow that A A = max„gN^. Ann — iiiin„gHj ^nn be 
much smaller than \Ajj\. Then, by the mean value theorem, the time limit (31) is again 
Ajj for any normalized Pnn(O). Since this is valid for any normalized distribution, it is 
possible to take, for simplicity, the uniform expression Pnn(O) = 1/^j, where Zj is the 
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number of states in N^. This expression corresponds to the microcanonical distribution, 
with the energy Ej. Since hmit (31) acquires the microcanonical form, one can say that 
there happens the cigcnstate thermahzation. For isolated quantum systems, the eigcnstate 
thermalization was suggested as a hypothesis by Deutsch [32] and Srednicki [33]. Here, 
we have shown that for a quasi-isolated system, with the prescribed initial condition, 
this is not a hypothesis but a rigorous asymptotic result. Note that the microcanonical 
distribution is a particular case of the general representative distribution [12,31]. 

It is possible, in principle, that, in addition to a measurable part of the density Pmn{^), 
there would exist a non-measurable term of the form 

Pmn(s) ^\Y1 ^rnnj[Sie - amnj) + S{e + amnj)] , 
3 

with positive coefficients Cmnj normalized so that the total density be normalized as in 
Eq. (25). This would yield the appearance in factor (27) of the fluctuating term 

-C'mn(^) = Cmnj COs{amnjt) ■ (36) 
3 

Then, instead of limit (31), we would have the asymptotic expression 

{A{t)) ~ Yl Pnn{0)Ann + J2 Prnn{t)AnmD^^{t) , (37) 

n m^n 

as time tends to infinity. In that case, the system does not equilibrate completely, but 
only partially, exhibiting permanent fluctuations. Though the temporal average of the 
latter expression is zero, the fluctuations can be rather strong and their existence may 
drastically change the system properties, hence these fluctuations cannot be neglected. 
A physical example of such a situation can be associated with heterophase fluctuations 
[29,34]. 



4 Irreversibility of time direction 

The problem of time irreversibility has attracted great attention (see, e.g., Refs. [21- 
23,35,36]). Usually, one connects the irreversibility with the behavior of the system ther- 
modynamic characteristics, such as the increase of entropy with time. The Gibbs entropy 
is known to be constant for isolated systems, because of which the Boltzmann entropy 
seems to be preferable for considering this problem [37,38]. 

Another possibility is to consider the average entropy or the average information. The 
latter is deflned as the average of the information content Inp^s, 

I{t) = TiABPABit) In Pas , (38) 

similarly to the deflnition of the operator averages (5). Then we have the following 
statement. 

Theorem 2. The average information (38) does not increase with time: 

I(t) - /(O) < . (39) 
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Proof. It is straightforward to see that 

T^ABPABit) In PAS = Tr ABPABU^it) {lnpAB)U{t) = 

^Tr^BPAfiln U'^{t)pABU{t)^ ^ TtabPab^t^ PAB{-t) ■ 

Therefore 

/(O) - I{t) = Tr AB [pab In pab - Pab In • 

Then we use the Gibbs-Klein inequahty, proved by Gibbs [39] for classical distributions 
and generalized by Klein [40] for quantum operators. According to this inequality, for 
any two non-negative operators A and B, acting on a Hilbert space T-L, one has 



Tr^ (A log i - i log > Trn{A - B) 



In our case, we get 
Since 



7(0) - I{t) > TiAB [pab - pAB{-t)] . (40) 



T^ABPAB{—t) — T^abPAB — 1 , 

we come to the sought inequality (39). 

The loss of information with time can be interpreted as the origin of the time irre- 
versibility. However, the irreversibility of time is a phenomenon common for all macro- 
scopic systems, whether thermodynamic or generally dynamic. Therefore the explanation 
of this phenomenon should not be based only on the notion of entropy. Another point 
of view on the time irreversibility, not necessarily connected with thermodynamics, relies 
on the fact that, strictly speaking, absolutely isolated systems do not exist, since any 
given system is always, at least weakly, connected to some surrounding [22,36]. A real 
system can be quasi-isolated, but cannot be isolated completely, since there always exist 
small perturbations, often uncontrollable, which induce the interaction of the system with 
its surrounding [41-43]. This point of view has been emphasized long ago by Borel [44]. 
Accepting this, we come to the conclusion that any given real system, at least partially, 
equihbrates, according to Eqs. (31) or (37). The fact that any real system tends to an 
equilibrated state implies that there exists irreversible evolution, that is, time is irre- 
versible. The partially equilibrated state can be a kind of a steady or quasi-steady state 
[45]. 



5 Conclusion 

The main results and conclusions of the present paper can be formulated as follows: 

(i) The temporal evolution of quasi-isolated quantum systems can be treated on the 
general level, valid for systems of arbitrary physical nature and under any initial condi- 
tions. 

(ii) Quasi-isolated systems, with increasing time, equilibrate, at least partially. 

(iii) When the system equilibrates to a stationary state, the latter, generally, keeps 
information on initial conditions and is characterized by a representative ensemble. 
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(iv) For some initial conditions, the eigenstate thermalization happens, being a rigor- 
ous asymptotic result. 

(v) Accepting that all real systems are at the most only quasi-isolated, but never 
absolutely isolated, defines the arrow of time, explaining why time is irreversible. 

The work is supported by a grant of the Russian Foundation for Basic Research. 
Discussions with E.P. Yukalova are appreciated. 
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